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ABSTRACT 

Polyampholytes (PAs) are polymers composed of (quenched) random sequences of 
positive and negatively charged monomeric groups. We show that the radius of gyration, 
R g , of a PA strongly depends on its overall excess charge Q, and is very weakly influenced by 
other aspects of the sequence. For Q < Q c « qoVN , where N is the number of monomers 
of charge ±qo, the PA is compact, while for Q > Q c it is stretched. Some aspects of this 
behavior can be understood by analogy with the shape instability of a charged drop where 
the Rayleigh charge plays the role of Q c . 
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As the building blocks of living matter, long chain macromolecules have been the 
subject of intense study. Whereas there is now a reasonably firm basis for understanding 
the physical properties of homopolymers considerably less is known about the heteropoly- 
mers of biological significance. From a biologist's perspective, it is the specific properties 
of a particular molecule that are of interest. After all, the genetic information is coded 
by very specific sequences of nucleic acids, which are in turn translated to the chain of 
amino acids forming a protein[l]. In response to the van der Waals, hydrogen bonding, 
hydrophobic/hydrophilic, and Coulomb interactions, the protein folds into a 'native' shape 
that is responsible for its activity. Given the large number of monomers making up such 
chains, finding the configuration of a particular molecule is a formidable task. By con- 
trast, a physicist's approach is to sacrifice the specificity, in the hope of gleaming some 
general information from simplified models [2]. There are in fact several statistical descrip- 
tions of proteins as random linear sequences of elements with a variety of interactions that 
determine their final shape [3] [4]. 

Statistical models only provide gross, general descriptions of the heteropolymer. Per- 
haps the coarsest characteristic of a molecule is some measure of its size, such as the radius 
of gyration, R g . For a polymer with charged groups, the overall size is most likely con- 
trolled by the Coulomb interactions which are the strongest and with the longest range. 
In this paper we consider the typical size of poly ampholytes (PAs): polymers composed of 
positive and negatively charged groups [5]. Even determining the size of a simple model 
PA, formed from charges ±qo randomly dispersed along the chain, has proved controver- 
sial. Appealing to a Debye-Hiickel (DH) attraction in analogy to electrolytes, Higgs and 
Joanny[6] conclude that a PA minimizes its (extensive) free energy by contracting to a 
compact structure. By contrast, in ref.[4], we used an analogy to the renormalization 
group (RG) treatment of homogeneously charged polymers [7] to conclude that the PA is 
stretched by the (non-extensive) energy cost of typical charge fluctuations. This contra- 
diction was partially resolved [8] by noting that the DH theory requires exact neutrality 
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of the electrolyte, while the RG-inspired approach treats equally all possible sequences; 
the typical charge of a random sequence, Q ~ ±qo\^N, sufficing to modify the behavior 
of the PA. Monte Carlo simulations [8] indeed confirm that PAs with Q = compactify at 
low temperatures, while sampling all random quenches with unrestricted Q produces an 
average R g oc N. Unfortunately, the broad nature of the distribution of sizes in the latter 
ensemble provides little information about R g of a specific PA. 

Recently, Yu et al.[9] performed a detailed investigation of the volume transition in 
gels produced by crosslinking PAs. The screening length in the experiments is quite big, 
enveloping a large number of monomers N. By changing the pH of the solution, it is 
possible to gradually modify the excess charge Q, within the screening length. The neutral 
gel is naturally the most compact. Surprisingly, the volume of the gel does not change 
appreciably in an interval of Q below a threshold Q c ~ qoN 1 ^ 2 . Immediately beyond 
the threshold, the volume of the gel increases by more than an order of magnitude. These 
results motivated us to make a systematic examination of the dependence of the R g of a PA 
on its excess charge. Interestingly, we find that Q is quite a good indicator of the overall size 
of a PA, i.e. PAs with quite different sequences, but the same excess charge, have a narrow 
distribution of R g s. The broad nature of the distribution of R g s for unrestricted polymers 
solely reflects the drastic variations of size with Q. In agreement with experiments, we 
find that PAs are compact for small Q, and stretched when Q exceeds a critical value of 
Q c ~ qo\^N. We provide a (partial) justification for this transition by analogy to a charged 
drop: The drop is stable at small Q, but disintegrates after suffering an instability when 
Q exceeds the Rayleigh charge Qr. 

The Monte Carlo (MC) procedure used in this work is identical to that of ref.[8]. 
We use self-avoiding chains on a cubic lattice of spacing a, with electrostatic interaction 
Uij{r) = %^/V2a 2 + r 2 for every pair of charges at a distance r, where = 

±go with signs specific to a particular quench. We examined the dependence of polymer 
size on both temperature, T, and excess charge, Q. Since we considered systems with 
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quenched randomness, the results were always averaged over 10 different quenches for 
each T and Q. The MC time unit is defined as the interval over which N attempts are 
made to move atoms. (Due to the long-range interaction, CPU time per MC time unit 
grows as N 2 .) Each quenched configuration at a given temperature was equilibrated for 
250iV 2 MC time units, which for N = 64 amounted to about 100 statistically independent 
configurations[10]. To collect all the data in Fig.l we needed about two months of CPU on 
a Silicon Graphics R4000 workstation. Equilibration times were sufficiently long to insure 
that statistical uncertainties of thermal averages were smaller than the differences between 
different quenches. 

Fig.l depicts the temperature dependence of R 2 for 64-monomer chains. The number 
near each curve indicates the charge, Q/qo- At very high temperatures the electrostatic 
interactions are unimportant and the chains behave as self-avoiding walks, with R g oc N v 
and v — 0.588. As T is lowered, the effects of interactions are first observed for the strongly 
charged polymers at Tq ps Q 2 /R g ~ Q 2 /N v (measured in energy units). This follows from 
the estimate of Q 2 /R g for the electrostatic energy of a typical configuration. However, 
for weakly charged polymers, with Q < qo\/N, the main interaction is between the non- 
homogeneities in the charge distribution. The energy of a typical density fluctuation, 
q(Q — q)/R g ps —q 2 N/R g for q ps qoVN exceeds Q 2 /R g , and such PAs start to deviate from 
their infinite-T behavior at Tq ps q 2 N/R g ~ q^N 1- ". These high temperature estimates 
suggest that on lowering T, the size of PAs with large Q increases (since the Coulomb 
energy of excess charge is repulsive), while R g decreases for PAs with small Q (since 
the charge fluctuation energy is attractive). This is confirmed by a (high temperature) 
perturbative expansion of the radius in the strength of the interaction as in ref.[8]. The 
first correction is proportional to q^qj = (Q 2 — q 2 N)/N 2 for i ^ j, and the change in the 
behavior of R g occurs for Q c = qoVN. (There is no a priori reason to expect higher order 
terms to also vanish at Q = Q c .) Are the averages in Fig.l a meaningful measure of the 
PA size? Fig. 2 depicts histograms of the distribution of R 2 at T = 0.2q 2 /a for several 
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values of Q. Since thermal fluctuations are small, this histogram represents the differences 
between quenches. Note that the distributions are fairly narrow; their widths not exceeding 
the distance between their averages. Thus the location of a point in Fig.l provides a 
good measure of R 2 (Q), without need for any additional reference to the details of the 
sequence. If this function is known, the average of R 2 for unrestricted quenches is simply 
obtained from, /^(random) = J °° dQR 2 g (Q)P(Q), where P(Q) oc exp[-Q 2 /(2^A)] is 
the probability density of an excess charge Q. In previous work [8], we found that the 
distribution of R 2 g (random) was very broad. Although the sample of quenches was small, 
it included several completely collapsed and strongly stretched configurations. The current 
results indicate that this follows simply from the strong dependence of R 2 on Q, rather 
than a large scatter of R 2 among different quenches with similar Q. 

The curves in Fig.l have a monotonic dependence on temperature. This suggests that 
the line separating compact and extended states in the {Q,T) plane is roughly straight, 
starting from Q = Q c = qoVN at infinite T. To test this hypothesis, we considered the 
Q and N dependence of the radius of gyration for chains of lengths N = 16, 32, 64, 128. 
To achieve good thermal averages, simulation were performed at T = 0Aq 2 /a and not at 
the lowest temperature in Fig.l. The dependence of R 2 on Q is depicted in Fig. 3; the 
vertical axis is scaled so as to remove the trivial A-dependence for Q = 0. The charges 
on the horizontal axis are scaled to Q/Q C {N) for all polymer lengths. Although R 2 is a 
monotonically increasing function of Q, its variations are not gradual: for small Q the 
radius barely depends on Q, while beyond a certain threshold an extremely fast increase 
begins. This feature is further amplified at low temperatures: the low charge polymers get 
smaller while the charged ones increase in size. Fig. 3 strongly suggests that the transition 
from compact to stretched configurations at low temperatures still occurs for Q m Q c . In 
fact, plotting R 2 {T = 0Aq 2 /a) / R 2 g (T = oo) as the function of Q/N 1 / 2 , we find that the 
curves become steeper with increasing N and intersect at Q/N 1 ^ 2 1.4go, implying that 
for Q > 1.4Q c the PAs at low T are more stretched than self-avoiding walks. 
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To explain the above results, we start with the empirical observation that the neutral 
PA compactifies to more or less a spherical shape. This suggests an energy (or more 
correctly a quench averaged free energy) of the neutral PA that grows with the number of 
monomers as E(Q = 0) = — e c N + AixR 2 ^, where e c oc q 2 /a is the energy gain per particle, 
and R xs aN 1 / 3 is the radius, in the collapsed state. A positive surface tension, 7 xs q 2 /a 3 , 
accounts for the spherical shape of such PAs. If we now uniformly add a charge Q to each 
configuration, its energy increases by Q 2 /R. Therefore, E(Q) < — e c N+4nR 2 'j+Q 2 / R; the 
inequality indicating that the system can reduce its energy by rearranging the charges. We 
would like to find the rearrangements that minimize the energy. For the sake of generality 
we shall consider Q oc qoN 13 , with emphasis on the most relevant case of (3 = \. To gain 
some insight into the behavior of PAs we shall explore analogies to charged drops. In 
the following paragraphs we briefly discuss the shape of a charged conducting drop. This 
analogy is most appropriate for an annealed version of the problem in which the charges 
are free to move along the polymer chain. We then go on to consider the shape of a 
charged insulating drop (immobile charges), and finally conclude with applicability of such 
arguments to quenched PAs. 

With the charge on the surface of a freely suspended conducting spherical drop of 
radius R, the non-extensive contribution to the energy is E'(Q) = AirR 2 ^ + Q 2 /(2R). 
As Q increases the shape of the drop may change. The surface energy of an uncharged 
drop, £"(0) = 47ri? 2 7, sets the overall energy scale of the problem, while the dimensionless 
parameter, 

a = Q 2 /(16^ 3 7) = Q 2 /Q 2 r , 

determines its shape. (Qr is the Rayleigh charge of the drop.) Note that the above 
estimates of 7 and R for the model PA lead to Qr ks Q c . We shall initially consider only 
deformations in shape that maintain the integrity of a single drop. Such deformations of a 
charged conducting drop, and related problems, have been considered by many authorsfll]. 
It was noted by Rayleigh[12] that for a > 1 the charged sphere is locally unstable and 
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must thus become distorted. The problem was also studied by Taylor [13] who considered 
distortions into prolate spheroids. For spheroids, the energy and shape can be determined 
by minimizing a trial energy with respect to the eccentricity, e. We can show[10] that 
for a > 0.899 the energy of the spheroid is smaller than that of the sphere, and the drop 
deforms into a strongly elongated shape. This qualitative similarity to the experimental 
results[9], and to Fig. 3 is encouraging. While the drop strongly elongates for finite a, 
it maintains a finite aspect ratio. Gutin and Shakhnovich[14] recently reached a similar 
conclusion by assuming that the charged PA distorts into a cylindrical shape with one long 
dimension R\\, and two short dimensions R± which are related by volume conservation: 
Minimizing estimates of E(Q) for such shapes leads to i?y ~ jv( 4/3_1 )/ 3 , i.e. no significant 
stretching for (3 = 1/2. Another recent study by Dobrynin and Rubinstein[15] relaxes 
the constant volume constraint and also reaches the conclusion that there is an onset of 
stretching for (3 = |, although a completely stretched state is reached only for f3 = 2/3, 
when the Coulomb energy becomes extensive. 

However, the spheroidal shape is not a local energy minimum[13] as the strongly 
elongated spheroid is unstable to a variety of perturbations [11]. Experimentally, for a > 1 
the drop disintegrates into smaller ones. Indeed if splitting of a continuum droplet is 
allowed, its energy can be reduced to E'(0). This is achieved by splitting away from the 
original drop an infinite number of infinitesimal droplets which contain all the charge but 
have vanishing total electrostatic energy, surface area and volume, and removing them to 
infinity. Even if the particles are required to stay connected (to maintain the connectivity 
of a chain), the system may be able to reduce its energy by expelling its charge in the 
form of a finger. Balancing the Coulomb energy (Q 2 /L) of a finger of length L, with a 
condensation energy e c L/a, indicates L oc Q. Such fingers appear spontaneously if their 
cost (roughly Q\/e c /a) is less than the Coulomb energy of the uniformly charged sphere, 
Q 2 /R, i.e. for Q > R^Je c /a. Thus the typical annealed PA will have a protruding finger 
of length L oc N 13 for f3 > 1/3. However, since the weight of the finger is small, it will not 
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effect the scaling of R 2 for f3 < 5/9. Such annealed PAs have large spanning sizes without 
appreciably greater R g . 

Now consider an opposite limit where the charge is uniformly distributed along the 
polymer, mimicked by a droplet of charged insulator. The sum of electrostatic and surface 
energies for a spherical drop is 3Q 2 /(5R) + AtkR 2 ^. Unlike the conducting case, energy 
cannot be lowered by expelling charge into small drops or fingers (since distortions of 
charge and volume are now connected). However, for a > 0.293 a drop can reduce its 
energy by splitting into two equal drops separated by an in infinite distance. Additional 
splittings occur for larger charges; the optimal number of droplets increasing linearly with 
a. We can again constrain the overall object to remain singly connected by attaching the 
droplets with narrow tubes of total length L and diameter a. As long as L <C aN, most of 
the charge remains in the spheres. Equating the total electrostatic energy, oc Q 2 / L, to the 
condensation cost, e c L/a (or a surface energy of •yaL in the continuum parlance), leads 
to L oc Q. Not surprisingly, the tubes in this case behave much the same as the fingers 
for the conducting polymer. However, whereas the surface tension in the conducting case 
results in one big central drop, for the insulating case the droplets are separated as in a 
necklace. The radius of gyration is now of the same order as the span of the necklace, L. 

How can we relate the above discussion to the quenched PA? The immobile charges are 
likely to produce a necklace of compact beads joined by narrow strings once the Rayleigh 
charge has been exceeded. However, it may be possible for the random PA to reduce 
its energy (compared to a homogeneous one) by taking advantage of charge fluctuations 
along its sequence: The strings will be constructed of segments with a larger excess charge 
than average, while the beads are from segments closer to neutral. Thus, R g oc Q may 
underestimate the size of the PA. The necklace picture may be hard to quantify if the beads 
and strings have a broad distribution of sizes. The following argument seems to suggest an 
upper bound to the total length of the strings L (and hence for R g ). The string segments 
must lose at least a portion of the condensation energy associated with compact beads. 
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Since this energy loss can only come from a reduction in the Coulomb energy, we have 
L < aQ 2 /(e c R) oc jV 2/3_1 / 3 . Thus while the quenched PA is certainly extended beyond 
the Rayleigh instability, it seems difficult to justify R g oc N. However, we began the 
discussion with rather robust assumptions on the condensation energy (oc N) and surface 
energy (oc iV 2 / 3 ), while the end-result involved energies of order iV 1 / 2 . Consequently, an 
accurate estimate of the stretching of the quenched PA requires a more detailed knowledge 
its energy spectrum. A full understanding of the behavior of PAs is yet to be attained. 

We would like to thank B.I. Halperin for bringing the importance of surface tension 
to our attention, and D. Ertas for helpful discussions. This work was supported by the 
US-Israel BSF grant No. 92-00026, by the NSF through No. DMR-87-19217 (at MIT's 
CMSE), DMR 91-15491 (at Harvard), and the PYI program (MK). 
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Figure Captions 

Fig. 1: Rg (in units of a 2 ) as a function of T (in units of q^/a) for several values of 
the excess charge Q for a 64-monomer chain. Each point is an average over 10 
quenches. Independent quenches are used at different temperatures. The numbers 
near each curve indicate Q/qo- 

Fig. 2: Histograms of the distribution of (ten) values of R 2 g (measured for N = 64 at 
T = 0.2(/Q/a) for several values of Q/qo, indicated near the histograms. 

Fig. 3: Scaled R 2 g as a function of scaled charge Q/qo for chain lengths N = 16 (open 
triangles), 32 (full triangles), 64 (open circles), and 128 (full circles). 
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